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INFRARED EMISSION FROM THE ATMOSPHERE ABOVE 200 KM 


Janet P . Simpson 
Ames Research Center 


SUMMARY 


Many of the atoms and molecules of the earth's atmosphere radiate in the 
infrared. This radiation is calculated from to lOOOp for altitudes between 
200 km and 400 km; only zenith lines of sight are considered. The excitation 
of the atoms and molecules is due to collisions with other molecules, and to 
absorption of radiation from the earth and sun. In some cases the abundances 
of the molecules had to be estimated. The most important lines are the for- 
bidden lines from atomic oxygen at 63. Ip and 147p, and the vibration-rotation 
band of nitric oxide at 5.3p. These lines can have intensities as high as a 
few times 10“^ ergs cm”^ sec~^ sr”^ at 200 km altitude. In addition, the 
vibration-rotation bands of NO"*" at 4.3p and CO at 4.7p and the pure rotation 
lines of NO and NO'*' could be detected by infrared telescopes in space. 


INTRODUCTION 


Infrared telescopes on space satellites will have to look through the 
tenuous gas of the earth's outer atmosphere. This gas itself radiates in the 
infrared, and the radiation from the gas will enter the detectors of the tele- 
scope along with the radiation from the stars, planets, and galaxies that are 
the goals of the observations. An attempt is made here to estimate the 
amounts of infrared radiation that will enter the shuttle- and satellite-based 
telescopes and their detectors from the overlying natural atmosphere. 

The gas that is left at this altitude is composed mainly of neutral 0, 

N 2 , O 2 , Ar, H, and He. Of these, H and He produce line radiation only through 
recombination from ionized H and He. The spectrum of Ar also contains no low 
excitation infrared lines. N 2 and O 2 , being homonuclear diatomic molecules, 
have no dipole moment and hence have no permitted rotational or vibration- 
rotational lines. However, O 2 has two unpaired electrons (ground state ^E) 
and hence has a forbidden magnetic dipole rotation spectrum. Neutral atomic 
oxygen has very strong forbidden transitions among the fine structure levels 
of its ^P ground term. In addition, there should be neutral species with 
abundances too small to detect; these species have diffused upwards from the 
lower atmosphere. There are no low excitation infrared lines of neutral N, 
but NO, CO, and CO 2 have strong vibration-rotation bands in the infrared. 

Their abundances will have to be estimated. 

Ultraviolet photons from the sun can ionize any of the neutral atoms and 
molecules just mentioned. Charge exchange reactions also affect the abundances 



of the ions. From the UV flux from the sun and the known abundance of neu- 
trals, it should be possible to calculate the abundance of ions and minor 
constituents. However, many reactions are involved and the cross sections 
and diffusion rates are often known poorly, if at all. Fortunately, measure- 
ments have been made of the abundances of the ionic species under close to 
average conditions. 


THE ATMOSPHERE ABOVE 200 KM 


The atmosphere above 200 km is thought to be in diffusive equilibrium 
for most constituents. Numerical models, such as those described in CIRA 1972 
(ref. 1), must be used to obtain the density of the atmosphere which is 
strongly dependent on the kinetic gas temperature at these altitudes. The 
main differences between the models in CIRA 1972 and those in the U.S. Stan- 
dard Atmosphere Supplements, 1966 (ref. 2) is that the former models contain 
more neutral oxygen and less O 2 and N 2 to conform better to observations at 
150 km. Consequently, the temperature required to produce a given density is 
also different. The tables in CIRA 1972 give the density of each component 
as a function of altitude for exospheric temperatures ranging from 500° K to 
2200° K. The temperature of the atmosphere increases very rapidly from 
210-270° K at 100 km to the final exospheric temperature. At 350 km the 
temperature is only about 3 percent less than the exospheric temperature. 

The models in the tables are static models; they do not take into account 
the dynamical variation of temperature and density of the atmosphere. The 
main temperature variations are due to the solar cycle, short-term changes in 
the active regions on the sun, and the daily variation due to the rotation of 
the earth, and the geomagnetic activity. Given the parameters for the solar 
and geomagnetic activity and the time of day and latitude on earth, CIRA 1972 
(ref. 1) gives formulas for computing the appropriate exospheric temperature. 

The mean exospheric temperature is found to be about 1000° K; the temperature 
is less than 550° K for only 1 percent of the time and greater than 1900° K 
for only 1 percent of the time. The daytime temperature is higher by a factor 
of 1.3 than the nighttime temperature for the same solar flux. Over a solar 
cycle, the number of days with each value of the exospheric temperature forms 
a distribution function that is double peaked if the daytime and nighttime 
temperatures are considered separately. The two peaks correspond to solar 
minimum and solar maximum. Nicolet (ref. 3) investigated the temperatures 
for the solar cycle that had its maximum in 1957-1958. For this cycle, the 
peaks of the distribution function occurred at about 900° K and 1500° K for 
daytime temperatures and 700° K and 1200° K for nighttime temperatures. It is 
now apparent that the peak temperatures of the last solar maximum were several 
hundred degrees lower and that the 1957-58 maximum was exceptionally high. 
Because of the day-night variation, a satellite will see the temperature vary 
by up to 30 percent (depending on angular distance to the point of temperature 
maximum) every half orbit and the density vary by up to a factor of two. 

(For diffusive equilibrium, the density N(,R) goes as exp(-i?/ff) where 
the scale height, H, equals kT/mg, k is the Boltzmann constant, T is the kinetic 
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gas temperature, m is the mass of the constituent, g is the acceleration of 
gravity, and R is the altitude. If the temperature is constant, one can inte- 
grate N(R) to get the column density = N^H exp[-(i? - Rq)/H] where is the 
density at Rq. Then the change in column density with temperature is 


A col, den. 
col. den. 



where Rq is the height at which the density is constant, about 90-120 km. Of 
course the temperature is not constant down to Rq = 120 km, but the above 
expression gives the correct dependences on i?, m, and T.) 

In addition to the density variations due to temperature, there is a 
semiannual density variation of unknown cause: the density can vary by up to 

15 percent from the densities in the tables for a given temperature at an 
altitude of 350 km. The maximum variation is over 50 percent at 800 km. 
Again, CIRA 1972 (ref. 1) gives formulas and tables for computing these 
variations . 


All the formulas for the temperature and density variations represent 
fits to the observed data. For any particular day they may be considerably 
off, but on the average they should be as accurate as the data. Since there 
is so much variation in the temperature and density of the upper atmosphere, 
it is impossible to follow the atmosphere through one cycle of a satellite's 
orbit until after the satellite has been flying, when the required parameters 
of the average solar flux, solar activity, and geomagnetic activity are known. 
In order to get representative values of the infrared radiation, three tem- 
peratures were chosen for this computation: the average 1000° K and the 

extremes 550° K and 2000° K. The densities corresponding to these tempera- 
tures were used to compute average and extreme radiation from neutral 0 and 
O 2 . The information on ions and neutral NO corresponds to the temperatures 
not far different from the average 1000° K. 


INFRARED RADIATION 


The specific intensity in a line or band I is given by 

I{R) = r°°dv ,(j) N expfe d B .N - S, dr) dr- 

Jq Jfi 4tt ur\) u ^\4tt Jp ul u lu I J 

where v is the frequency, h is the Planck constant, and Bj^ are 

Einstein coefficients, ^y/^total ^Z/^total the fractions of atoms or 

molecules in the upper and lower levels respectively, and is the normalized 
line profile function, ^total the number density of a particular species 
of atom or molecule. Above 200 km, the lines are optically thin and I is 
approximated by 

hv 

^ A dr 
, 4ir ul u 
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which becomes 


I(R) 


hv 


A In 


4-rr ul u} total 


(col. den.) 


if the excitation is constant with altitude R. Col. den. is the column den- 


sity ^total atoms or molecules/cm^ . 


The intensity I will usually be expressed in ergs cm~^ sec“^ sr~^. It 
is converted to W cm~^ sr“^ by multiplying by 10“^. Wavelengths are given 
in microns, y. The results of the calculations for each molecule and line 
will be plotted in figures l(a)-(e). 


Neutral Atomic Oxygen 

The ground state of neutral oxygen is a triplet of the 2p^ configura- 
tion. Fine structure transitions within the ground state produce forbidden 
lines in the far infrared at 63. Ip (%2 ~ a^id 147p (^Pj - ^Pq) • The 

familiar optical forbidden lines of neutral oxygen are at 5577 A (^^2 ~ 
and 6300 A (^?2 - ^P2) • The upper levels are excited and deexcited by colli- 
sions and radiation. The equilibrium populations are given by the Boltzmann 
equation N 2 IN 1 = QzIQi exp[-(P 2 “ where is defined by this equa- 

tion and gi and Q 2 are the statistical weights of levels one and two, respec- 
tively. If Pgjj is the same as the translational temperature of the neutral 
gas, T, the populations are in local thermodynamic equilibrium (LTE) . Saraph 
(ref. 4) has given cross sections for collisional excitation by electrons. 
Calculations with the observed electron densities show that the level popula- 
tions are within 90 percent of the LTE population at 350 km and within 80 per- 
cent at 500 km for T = 1000° K for excitation by electrons alone. Collisions 
with neutral atoms and molecules will increase the populations of the upper 
levels. In fact, the populations of the upper levels could be greater than 
the LTE populations, because the electron temperature is substantially higher 
than the gas temperature. (At 550° K, the populations may be less than the 
LTE populations because the gas density is so small) . For the calculations 
below, it was assumed that the upper levels are populated according to the 
Boltzmann equation at the temperature of the neutral gas because the density 
of the neutrals is two to three orders of magnitude greater than the electron 
density. In any event, the LTE populations do not vary greatly for the tem- 
peratures under consideration (see table 1). 


TABLE 1.- LTE POPULATIONS FOR THE GROUND STATE OF NEUTRAL OXYGEN 


T, °K 


17(3Pi)//7^^t 


550 

0.073 

0.263 

0.664 

1000 

.089 

.294 

.617 

2000 

. 100 

.314 

.586 

3000 

. 101 

.335 

.563 
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For the calculation of the [01] line intensities, the number densities 
were taken from CIRA 1972 (ref. 1) and the Einstein coefficients from Garstang 
(ref. 5): A^Pi - = 1.7x10-5 sec“^ and /1(3P2 ~ = 9.0xl0"5 sec"!. 

The results are given in table 2. The optical depth in the line center 
t(Av = 0) is given by 


A 


O/5)N0p^) - il/(3Pj) 


63.1y(Av-0) 8-n(2vkT/m)^^ ^ 


X col. den. 


if there is no turbulent line broadening. These optical depths are so small 
at 300 km that self-absorption is negligible. However, at low altitudes the 
optical depth is large enough that self-absorption should be considered 
(Kockarts and Peetermanns, ref. 6). 

The rocket observations of [Ol] that have been made have all been for 
altitudes below 120 km (Feldman ^nd McNutt, ref. 7). 


TABLE 2.- LINE INTENSITIES FOR NEUTRAL OXYGEN 


Altitude , 

Column density. 

P, ergs cm ^ sec ^ sr ^ 


km 

-"2 

cm 

63. Ip 147p 

^ Q 9 U 3 • 1 ^ 

T = 550° K 




1.52x1033 

6,86x1033 

3.22x103'' 

1.58x1035 

8.29x1035 

9.0x10 

4.1x10 

1.9x10 

9.4x10 

4.9x10 


= 1000° K 


= 2000° K 




9. 18x1035 

6.5x10 

1.39x1035 

9.8x10“'' 

2. 14x1035 

1.5x10“3 

3.39x1035 

2.4x10-3 

5.79x1036 

4.0x10-3 




2.0x10-3 

0.00001 

9.2x10-3 

.00006 

4.3x10~7 

.0003 

2.1x10~3 

.001 

1.1x10-5 

.008 


7.97x1q3'' 

5.3x10-5 

1.3x10-6 

1.83x1035 

1.2x10"'' 

3.0x10"6 

4.30x1035 

2.9x10"'' 

7.0x10"6 

1.05x1036 

6.9x10 *' 

1.7x10-5 

2.72x1036 

00 

X 

o 

CO 

4.3x10-5 


1, 

,7x10"5 

CM 

,5x10"5 

3. 

,9x10-5 

6, 

. Ixl0"5 

1. 

,0x10"'' 


















©2 and N 2 


The homonuclear diatomic molecule, O 2 , has no electric dipole moment. 

Thus the electric dipole vibration-rotation and rotational spectra are 
strictly forbidden. However, there are two unpaired electrons in O 2 , giving 
a ground state of and a magnetic dipole moment. The resulting magnetic 
dipole rotational spectrum has been observed and discussed by Gebbie et at. 

(ref. 8). These lines have also been observed in the solar absorption spec- 
trum. Gebbie et at. give theoretical and observed values for the integrated 
line strengths iS = /k dv where ic = dx/dr". The integrated line strength is 
proportional to the Einstein A coefficient times the population of the lower 
energy level. For the strongest lines between 165y and 705y, the maximum 
theoretical values of A are of the order of 10“^ sec~^; the observed values 
agree within 25 percent. The observed values of A were used and ^4 = 10 sec ^ 
was assumed for transitions between the unobserved higher rotational levels. 

The molecule O 2 is described by Hund's case b; the observations of Gebbie et at. 
only cover wavelengths with transitions as high as quantum change t\K = 11 - 9, 
but rotational energy levels were calculated to K = 27. The population of a 
level is given by 


N(J,K) ^ (2J + 1) exp[-7i:(js: + l)hoB/kT] 

^ -hkTllhoB 
tot 

if the energy levels are functions of K only (actually there is a small split- 
ting with t7) . The rotational constant is B (McKnight and Gordy, ref. 9) and 
J = 7([-l, K, ox K 1. For T = 300'’ K, the maximum population occurs at 
K = 7 , J = 8, but for T = 1000'’ K the maximum population ±s at K = 15 - 17, 
and for T = 2000° K for the maximum is at = 21 - 23. (Only odd numbered 
values of K are populated.) Therefore, at 1000° K estimates of the inten- 
sities for the strongest lines require extrapolation to get values for the 
Einstein A coefficients. 

In table 3 the transition = 23 - 21 is the strongest for T = 1000° K 
but = 27 - 25 is strongest for T = 2000° K. Normally the Einstein A coeffi- 
cients increase with decreasing wavelength; however the O 2 rotational spectrum 
becomes increasingly forbidden with increasing K (Gebbie et at. , ref. 8). 

The above estimates for the Einstein A coefficients and the column den- 
sities from CIRA 1972 (ref. 1) were used to predict the intensities of the 
strongest O 2 rotational lines in table 4. These lines are much weaker than 
the lines from neutral oxygen. 

The molecule N 2 has a much larger column density than O 2 (see table 4) . 
However, since N 2 is a homonuclear molecule with no magnetic dipole moment 
due to electrons, its rotation and rotation-vibration spectrum are strongly 
forbidden. Because the nuclear spin of the nitrogen atom is one, there is a 
nuclear magnetic moment for N 2 . The transitions, if they existed, would have 
an extremely small transition probability because of the small nuclear mag- 
netic moment; the intensities of any N 2 rotational lines would be much smaller 
than the intensities of O 2 rotational lines. 
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TABLE 3.- SOME ROTATIONAL LINES FOR O 2 


hK, AJ 
Transition 

X 

y 

A, sec ^ 

ho A J - 1 _ 1 

4irX N ’ sec sr 

tot 

T = 1000° K 

T = 2000° K 

3-1, 2-2 

705.8 

5.4x10"® 

7.4x10“26 

2.9x10“26 

5-3, 4-4 

387.4 

6.3x10"® 

2.7x10-25 

1.0x10“25 

9-7, 8-8 

204.4 

1.2x10"^ 

1.6x10“2'+ 

6.5x10"25 

11-9, 10-10 

165.4 

9.2x10“® 

1. 7x10“2‘+ 

7.6x10“25 

15-13, 14-14 

120.5 

1.0x10“^ 

2.9x10“2'+ 

1.4xl0"2't 

19-17, 18-18 

94.5 



3.4x10“2'+ 

1.9x10"2'+ 

23-21, 22-22 

77.8 



3.6x10“2‘+ 

2.4x10-24 

27-25, 24-24 

66.1 



2.9x10-2'+ 

2.7x10-24 


TABLE 4.- INTENSITIES FOR ROTATIONAL LINES OF O 2 


Altitude, 

km 

N 2 Column density, 
cm“2 

O 2 Column density, 
cm" 2 

Intensity, 

ergs cm"2 sec"l sr“l 

77. 8y 

66 . ly 

T = 1000° K 

400 

1.45x1013 

4.19x1011 

1.5x10-12 

1.4x10-12 

350 

6.23x1013 

2.21x1012 

7.9x10-12 

7.1x10-12 

300 

2.81x1014 

1.24x1013 

4.4x10-11 

4.0x10-11 

250 

1.35xl0l® 

7.67x1013 

2.8x20-1® 

2.5x10-1® 

200 

7.31x101® 

5.16x1014 

1.9x10-3 

1.7x10-3 

T = 2000° K 

400 

l.Olxlol® 

5.28x1013 

1.3xl0“l® 

1.4xl0"l® 

350 

2. lOxlol® 

1.23x1014 

3.0xl0"l® 

3.3xl0"l® 

300 

4.52x1015 

2.94x1014 

7.2xlO"10 

8 . 0 x 10 "!® 

250 

1.03xl0l® 

7.52x1014 

1.8x10-3 

2.0x10-3 

200 

2.66x101® 

2.24xl0l® 

5.5x10-3 



6.1x10-3 


Nitric Oxide 

Nitric oxide has not been observed in the earth's atmosphere at heights 
greater than about 150 km; calculations of the emission from NO must employ 
theoretical models and extrapolation. Ghosh (ref. 10) made a model for the 
NO abundance to 280 km; however, according to Strobel (ref. 11), Ghosh's 
reaction rates are in error. Strobel, Hunten, and McElroy (ref. 12) and 
Strobel (ref. 11) have also computed models for the NO abundance up to 
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240 km. The models of Strobel et al. are for an exospheric temperature of 
750° K and a variety of assumptions about the reaction rates. Strobel 
(ref. 11) computed the daily variation for a minimum exospheric temperature of 
940° K (that is, the maximum temperature in his model is -1200° K) . 

Column densities were estimated from the models of Strobel et at. 

(ref. 12) and Strobel (ref. 11). The density of NO was first extrapolated 
from 240 km to infinity. Above 120 km, the number density of a molecule is 
generally determined by diffusive mixing, although chemical reactions are 
important for NO lower than about 250 km. For diffusive mixing. 


NiR) = exp[-(i? - Ro)lE] 

where is the density at Rq. Obviously, the scale height ff(= kT/mg) is not 
constant with altitude. However, if E is constant, the equation for N can be 
integrated directly to get 

Col. den. (7?) = RJl exp[-(i? - Rq)/E] 


Values of 77 at i? = 300 km and 250 km are given in table 5. 

TABLE 5.- SCALE HEIGHTS FOR DIFFUSIVE MIXING 


Molecule 

Scale height , km 

T = 1000° K 

T = 750° K 


O 2 

29.26 

21.95 

i? = 300 

NO 

31.21 

23.41 


1n2 

33.44 

25.08 


f02 

28.76 

21.57 

7? = 250 • 

NO 

30.68 

23.01 


[n2 

32.87 

24.61 


However, since the temperature, and thus the scale height, are variable with 
altitude, it is better to use the scale heights for N 2 and O 2 for a complete 
model, and take the average to get the scale height for NO (table 6). Since 


TABLE 6.- SCALE HEIGHTS FOR NO 


T, °K 

= 230 km^, cm ^ 

E , km 

750 

1.6x10^ 

22.50 

940 

2.0x10® 

28.82 

1200 

6.0x10® 

36.79 
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Strobel (ref. 11) used models that are almost identical to the model in the 
U.S. Standard Atmosphere Supplements 1966 (ref, 2), scale heights from these 
models were used. Values of E and Nq were interpolated for T = 1000° K and the 
column density was computed. The column density for T = 1000° K is larger by 
a factor of 2 than the column density for T = 940° K at 350 km altitude. 

The errors in the column densities may be very large because of uncer- 
tainties in the models and the method of extrapolation. The model by Ghosh 
(ref. 10) is an order of magnitude less dense in NO at 230 km than Strobel' s 
(ref. 11) model. The column density at 85 km has been measured by Barth 
(ref. 13) to be 1.7xl0^‘* cm~^. However, the column density for Strobel's 
(ref. 11) model at 85 km is greater than 6x10^** cm"^. Strobel compared his 
results with the observations of Meira (ref. 14) between 70 km and 115 km. 

The integrated column density is larger by a factor of 2 for Strobel's 
calculations, but the number densities at 100 km (10® cm”®) agree "remarkably" 
well. Strobel noted that if the scale height for NO above 100 km were lower, 
the column densities would agree better, but then the density and column 
density at 250 km would be greatly lower. Strobel et at. (ref. 12) also com- 
pared their results to various observations from 85 km to 150 km. They con- 
cluded that the agreement is good. Rusch (ref. 15) also measured the density 
of NO to be ~10® cm ® at 105 km. Although atmospheric models vary greatly in 
density and temperature at high altitudes, such variation is very small around 
100 km. Thus the discrepancy between the observed and theoretical NO densi- 
ties cannot be ascribed to atmospheric variation. 

The vibration-rotation bands lie at 5.3y for the fundamental (Ay = 1), 

2.7y for the first overtone (Ay = 2), and 1.8y for the second overtone 
(Ay = 3), The wavelength can be considered to be constant for the band. The 
Einstein A coefficients are not measured directly; instead one measures the 
integrated absorption coefficient S = J k dv/Pc, where P is the pressure in 
atmospheres. For NO at STP, 0.999876 of the molecules are in the ground 
vibrational state and 


Ay=i 


^01^0 “ ^10-^0 and 


^ 02^0 °" ^ 20^0 


that ±s, A = 2.8038x10”® S/X^ cm“^ atm”^. The theoretical calculations of 
Michels (ref. 16) give Sqi = 134 cm“^ atm”^ (-^lo ” 13.2 sec”^); 

Sq 2 = 2.0 cm“2 atm”^ (^20 ~ 0.78 sec”^); and Sq^ = 0.08 cm”^ atm”^ 

(430 = 0.07 sec”^). Measurements of Sqi range from 74 cm”^ atm”^ to 
138 cm”^ atm”^ (ref. 17) with an average of 125 to 127 cm”2 atm”^. Measure- 
ments of Sq 2 range from 2.0 to 2,8 cm"^ atm"^. Since A S, one can quickly 
adjust 4 if a different value of S is preferred. 


The next problem is choosing the proper temperature to use in calculating 
the level populations. The rotational energies come into equilibrium very 
rapidly with the kinetic energies of the most abundant species, which are the 
neutral atoms and molecules. However, the cross sections for excitation of 
the vibrational levels are much smaller than the cross sections for excitation 
of rotational levels and one must compare in detail the excitation rates with 
the radiative depopulation rates. 
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A vibrational level of a molecule can be excited by radiative absorption, 
by radiative cascade, or by collisional excitation. The radiative absorption 
rate is given by f dv where is the mean intensity. Since 

^Zu ~ ^ vibration-rotation band, the radiative absorption 

rate due to the earth, is given by 

-i 

exp(/zv/fer^) - 1 


is the Planck function at the temperature of the earth, T^. The radiative 
absorption rate for the sun, Rq, is given by 


ul 


R^(t,u) exp(hv/kT ) 


- 1 


where = 6.800x10 ^ ster and Tq ~ 5850° K at 5p (ref. 18). Radiative cas- 
cade occurs when a higher vibrational level is excited (by radiation from the 
sun) and the level decays by cascade instead of returning directly to the 
ground state. This is important for vibrational levels because 
>> ^Ay=2 For level y = 1 the excitation rate from radiative cas- 

cade through u = 2 is given by 


A 


R. 


21 


i21 


+ A 


20 


[i?^(0,2) + i?^(0,2)] = i?^(0,2) 


There are two kinds of collisional excitation: (1) translation-vibration 

(,T-V) , where the kinetic energy of the colliding molecule goes into vibrational 
energy; and (2) vibration-vibration (V-V) where the vibrational energy of the 
colliding molecule is transferred to the other molecule. The second reaction 
is strongly resonant for molecules with similar vibrational energy levels. 

The reaction rates are usually expressed for collisional de-excitation for 
T-V and in the exothermic direction for V-V; that is, 


R, {T-V) = k, {T)NN 
deex deex NO coll 


and 


R 


ex 


iV-V) = 


k ^coll^" 

exo NO — 


or 


^coll(total) 


%(y=l) 

R. (V-V) = k N 

NO(total) 


coll 


where /l^coll number density of the colliding molecule (or atom, for T-V ) . 

For excitation, the rate coefficients k are related by exp(-/zv/kT) 
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and ~ ^exo exp (~?zv/?cT) . The vibrational levels are deexcited by spontaneous 

emission, by stimulated emission, by absorption to higher levels, and by col- 
lisions. However, in all the cases considered here, the spontaneous emission 
rate is much larger than the stimulated emission, absorption, or collisional 
rates, which thus can be ignored. Then the population of the first vibra- 
tional level is given by 



where Tgx is the excitation temperature. (Similar equations hold for the 
higher vibrational levels. The intensity in the band is equal to the sum of 
the intensities from each vibration rotation transition with Au = 1 . However, 
the 1-0 transition is by far the strongest.) The resulting excitation tem- 
perature of NO ranges from 261° K at 400 km to 284° K at 200 km on the shaded 
side of the earth, and from 272° K to 289° K on the sunlit side of the earth 
for T = 1000° K. For T = 2000° K, the excitation temperature is only 4° K 
higher. The small variation in temperature shows that radiation from the 
earth dominates the excitation process. Thus the excitation temperature would 
be higher or lower if the temperature of the surface of the earth were higher 
or lower than average. Estimated vibration-rotation band intensities are 
given in table 7 for T = 1000° K. Similar calculations of the excitation of 
NO have been done by Bauer et at. (ref. 19), Degges (ref. 20), and Kummler 
and Bortner (ref. 21). 


TABLE 7.- VIBRATION-ROTATION BANDS OF NO FOR T = 1000° K 


Altitude , 
km 

Column density, 
cm ^ 

■^5.3p. ergs cm 

2 - 1 - 1 

^ sec sr ^ 

Night 

Day 

400 

3.5x10^° 

4.4x10-2 

6.7x10-2 

350 

l.SxlO^l 

2.4x10-3 

3.5x10-3 

300 

9.3xl0ll 

1.3x10-3 

1.9x10“3 

250 

4 . 7 x 10^2 

8.3x10-3 

l.lxlO-'* 

200 

2.4x10^3 

6.5x10“'+ 

7.8x10“'+ 


The vibration rotation bands are spread out over a finite wavelength 
Interval that increases with increasing rotational temperature. Carpenter and 
Franzosa (ref. 22) have calculated the absorption spectrum for the y = 0-1 
transition of NO for temperatures of 200° K to 4000° K. If the individual 
rotation lines are unsaturated, the absorption widths for the whole band 
should be similar to the band widths in emission. The half-power band widths 
and the tenth-power band widths have been estimated from the graphs of Carpen- 
ter and Franzosa (ref. 22) and are given in table 8 as AA/A. The energy is 
quite spread out; narrow band spectrophotometry may not be bothered by the NO 
emission. 
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TABLE 8.- BAND WIDTHS FOR NO 



T = 1000° K 

T = 2000° K 


^^(1/2) 



A 


Fundamental 

0.10 

0.13 

First overtone 

.05 

.06 


> 

0 



A 


Fundamental 

0.14 

0.19 

First overtone 

.07 

.09 


Since NO has a permanent dipole moment, It also has a rotational line 
spectrum in the far infrared. Hall and Dowling (ref. 23) have measured the 
wavelengths of the lines and the rotational constants of the molecule. The 
ground electronic state is the doublet ^^^^2 the 

lie about 121 cm~ ^ above the ^^\/2 levels. The rotational spectrum consists 
of transitions with hJ = \ within the ^Hj^2 level and with AJ = 1 within 
the ^ng ^2 level. The number of molecules in an energy level J is given by 




J 


N 


tot 


= -^ (2J + l)exp 



where the partition function Q equals kT/hcB if the energy levels are given by 
Ej = BJ (J + 1) and the effects of higher order vibrational levels on the 
energies Ej have been neglected. Q\ partitions the molecules into 
and ^ 113/2 

ei(^ni/2) 1 + exp(-;zeEj2i/^^) 


exp {-hoE -^21 

^ 3 / 2 ^ 1 + exp {-hcE -^21^^^^ 

where 5’i21 1® the energy separating the ^^ 1/2 ^^3/2 levels (121 cm“^). 

The value of the transition probabilities j were estimated from the 

dipole moment. Townes and Schawlow (ref. 24) give the rotational transition 
probability as 

A — 64it** 2 =1+1 

^J+\,J ~ 2J + 3 
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where A is the wavelength of the line and y is the dipole moment of the mole- 
cule. Townes and Schawlow give the value for the dipole moment y = 0.16 debyes 
for J = 0. However, the dipole moment is not constant with J, and particularly 
will vary if the higher order terms in the energy level equation 
Ej = BJ {J + \) - DJ^{J + 1)^ + (J + 1)^ + . . . are important. Since the 
higher order terms are important for NO, it is thus only an approximation to 
use y = 0.16 debyes for all values of Aj^i j. 

The specific intensities of the NO rotational lines were calculated using 
the column densities for NO at 1000° K from table 8 and are given in table 9. 
Only a few lines are tabulated here since many of the lines have similar 
intensities because A/\ increases as decreases. (At T = 1000° K, the 

maximum level population occurs at J = 27/2, but the maximum emission occurs 
att/~61/2.) At 1000° K the intensities are small, but they should be of 
the order of 10^ larger for T = 2000° K. 


TABLE 9.- ROTATIONAL LINES OF NO AT 1000° K 


Altitude, 

km 

Intensity, ergs cm ^ sec ® 

sr“l 

109y 

J = 55/2 - 53/2] 

107y [^Il 3 / 2 > 

J = 55/2 - 53/2] 

82. 4y [^ni/ 2 » 

J = 12>H - 71/2] 

400 

1.8x10“® 

1.5x10“® 

1.6x10“® 

350 

9.2x10“® 

7.5x10“® 

8.3x10“® 

300 

4.7x10“® 

3.8x10“® 

4.3x10“® 

250 

1.0x10-7 

8.5x10“® 

9.4x10“® 

200 

8.7x10-7 

7.9x10“7 

7.4x10“7 


CO and CO 2 

Theoretical models of the abundances of CO and CO 2 have been calculated 
up to 200 km by Hays and Olivero (ref. 25) and up to 160 km by Hunt (ref. 26). 
The abundances below these altitudes are determined by diffusion and by the 
dissociation of CO 2 to CO. Since CO 2 is such a heavy molecule, its abundance 
falls off more rapidly than CO with increasing altitude. At 110 km both 
models give the density of CO 2 between 3x10® and 5x10® molecules/cm® and 
CO an order of magnitude less, but at about 160 km, the abundances of CO and 
CO 2 are equal. Densities and column densities were estimated for higher 
altitudes under the assumption of pure diffusive mixing. The molecule CO is 
very stable against ionization or dissociation, and diffusion should be 
accurate at high altitudes. Since CO has the same mass as N 2 , under diffusion 
it should scale as N 2 . The ratio of the CO abundance from Hunt's model 
(ref. 25) to the abundance of No for T = 1000° K increases from 2.6x10“® at 

— r ^ 

110 km to 8.5x10 ^ at 160 km because diffusion is not yet dominant at these 
altitudes. However, above these altitudes eventually diffusion dominates, 
and it was assumed that the column density of CO is equal to 8.5x10“® times 
the column density of N 2 (table 4) . 
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The vibration- rotation spectrum of CO consists of a fundamental band at 
4.67y and overtones at 2.35y and shorter wavelengths. The collisional 
excitation rates for CO and CO 2 have been intensively studied because of the 
use of the molecules in lasers; the values used here came from the review by 
Taylor (ref. 27). At T” = 1000° K, collisional excitation is relatively unim- 
portant compared to radiation, and the excitation temperature of the upper 
vibrational level of the 4.67y transition of CO ranges from 291° K at 400 km 
to 297° K at 200 km for molecules in the sun and from 264° K to 277° K for 
molecules shaded from the sun. However, at f = 2000° K the excitation tem- 
perature of CO ranges from 320° K at 400 km to 395° K at 200 km in the sun 
and from 311° K to 393° K in the shade. The intensities of the 4.67y 
vibration-rotation fundamental of CO are given in table 10 for T = 1000° K. 
For higher temperatures, the intensities should be increased proportionally 
to the increased column density of N 2 times the increased population from the 
higher excitation temperature. 


TABLE 10.- INTENSITIES FROM CO AT T = 1000° K 


Altitude, 

km 

Intensity, ergs cm ^ sec ^ sr ^ 

4.67y (night) 

4.67y (day) 

400 

1.2x10-8 

3.5x10-8 

350 

5.1x10-8 

1.5x10-7 

300 

2. 5x10- ■7 

6.8x10-7 

250 

1.3x10"8 

3.5x10-8 

200 

9.7x10-8 

2.1x10-8 


The pure rotational lines of CO are fairly weak because of the small 
dipole moment (0.112 debyes) . At T = 1000° K, the strongest line is at 87y. 

Its intensity is about 3.8xl0~^®ergs cm~^ sec”^ sr kt T = 2000° K, the 
strongest line is at 63y and has an intensity of about 

l.GxlO"® ergs cm~^ sec~^ sr“^. The total intensity of the CO rotation lines 
is almost two orders of magnitude less than the intensity of the NO rotation 
lines . 

CO 2 has no pure rotational lines but it has many vibration rotation bands, 
the strongest of which are the V 2 fundamental at 15y and the fundamental 
at 4.26y. At T = 1000° K, collisional excitation of the 15y band is negli- 
gible (^ex “ 198° K) but at T = 2000° K the excitation temperature increases 
to 232° K at 200 km. Radiative excitation from the sun is negligible. How- 
ever, at 4.26y, radiative excitation from the sun is most important. In the 
sunlight, the excitation temperature ranges from 297° K at 400 km to 301° K 
at 200 km; in the shade, however, the excitation temperature ranges from 228° K 
at 400 km to 260° K at 200 km. The intensities, given in table 11, are still 
very small because of the small column densities. 
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TABLE 11.- INTENSITIES FROM CO 2 AT T = 1000° K 


Altitude, 

km 

Column density, 

-2 

cm 

Intensity, 

ergs cm 2 sec 1 

sr-l 

4.26p (night) 

4.26y (day) 

15y 

400 

3.4x102 

1.2x10“15 

4.3xl0"l‘* 

d.dxio-i** 

350 

3.7x103 

i.dxio-i** 

CO 

1 

0 

X 

9.6x10-13 

300 

4.0X10"* 

1.7x10-13 

5.1x10-12 

1.0x10-11 

250 

4.4x10^ 

3.0x10-12 

5.7x10-11 

1.1x10-1° 

200 

4.8x10^ 

1. ixio-10 

0 

1 

0 

X 

0 

1.2x10-° 


Ions (0+, N+, 0++, N 2 '*‘, 02 "^, NO+) 

Ionic densities are determined by diffusive equilibrium only at high 
altitudes, and are determined by photoionization and chemical reactions with 
ions and neutrals at lower altitudes. (The transition region is between 
300 km and 500 km, depending on the density.) The ionic densities have been 
measured by ion mass spectrometers on various rocket and satellite flights. 
Unfortunately, the published densities above 300 km were all measured during 
quiet sun years, when T is less than average T = 1000° K and the emitted 
intensities are less than average. Holmes et at. (ref. 28) measured the ion 
densities on a day flight {T ~ 750° K) and a night flight (7 ~ 650° K) in 
1963. Brinton et dl. (ref. 29) measured the ionic densities on a day flight 
in 1966 (? ~ 850° K) . Hoffman et dl. (ref. 30) also measured the ion den- 
sities on a day flight in 1966 {T = 940° K) . This rocket flight coincided 
with the pass of the Explorer 31 satellite at 970 km, on which there were 
experiments that measured the kinetic temperature of the ions as 2100° K and 
the electron temperature as 3650° K and 3700° K. The ionic column densities 
from the measurements with the largest gas temperature (Hoffman et dl. , 
ref. 30) are given in table 12. 


TABLE 12.- OBSERVED IONIC COLUMN DENSITIES AT T = 940° K 


Altitude, 

km 


N+ 

ions/ cm2 

NO+ 

ions/ cm2 

Na'*' 

ions/ cm2 

1 

O 2 + 

ions / cm2 

0++ 

ions/ cm2 

400 

3.2x1012 


1.1x101° 

1. Ixl0l° 

8.5x10° 

5.5x10° 

350 

4.2x1012 

1.7xl0ll 

■TTrini*! 

2. lxlol° 

1.3xl0l° 

5.9x10° 

300 

5.5x1012 

1.9xl0l 1 

llfflnM 

4.2xl0l° 

2.5xl0l° 

6.3x10° 

250 

7.3x1012 


2 . 9 x 1011 

8.4xl0l° 

1. Ixloll 

6.7x10° 

200 

8. 3x1012 

2. Ixioll 

1.0x1012 

1.5xl0ll 

4.6x1011 

7.0x10° 


Of the above ions, 02 "*" and N 2 ''’ are homonuclear and have no rotational 
lines or vibration rotation bands in the infrared. The ion O"*" has a line at 
505y, but it is a transition between levels that are 3.3 eV above the ground 
state. With this small level population and its small transition probability 
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Garstang, ref. 5), the line should be no stronger than 3x10 ergs cm ^ sec ^ 
at 250 km, which is very weak indeed. 

The ions IT*' and 0"*^ are 2p^ configurations with ground terms. The 
lines of are at 203. Ty (^Pg - ^Pj) and 121. 7ii (^Pj - and the lines of 

are at 88.18y (^Pq - ^P\) and Sl.lly (^P^ - ^P£) • The level populations 
are given by the Boltzmann equation with the temperature equal to the temper- 
ature of the electrons since the cross sections for collisional excitation 
of 0++ and IT*" by electrons are sufficiently large. With the use of the transition 
probabilities from Garstang (ref. 5) and the column densities from table 9, 
the intensities for T = 940° K are given in table 13. 


TABLE 13.- INTENSITIES FOR N+ AND 0++ AT P = 940° K 


Altitude , 
km 


Intensity, ergs 

-2 - 1 - 1 
cm ^ sec ^ sr ^ 


[N II] 121. 7y 

[N II] 203. 7y 

[0 III] 51.71y 

[0 III] 88.18y 

400 

0 

1 — H 
1 

0 

X 

7.9x10“^^ 

8.9x10-10 

8.9x10-11 

350 

0 

1 

0 

X 

0 

o> 

9.3x10-1^ 

9.5x10-10 

9.5x10-11 

300 

1.0x10-^ 

i.oxio-io 

l.OxlO-0 

1.0x10-10 

250 

1. 1x10“^ 

1. Ixl0"10 

1. lxlO-0 

1. lxlO-10 

200 

1. Ixl0~9 

1.2x10-10 

1. lxlO-0 

0 

1 — 1 

1 

0 

X 


The line strengths for NO"*" have not been measured in the laboratory 
because NO"*" is so reactive. However, Billingsley (ref. 31) has calculated 
Einstein coefficients for the vibration-rotation fundamental (4.27y) and first 
overtone (2.15p) and the dipole moment as a function of vibrational level. 

Line NO, NO^ is inf rared- active and radiates in its vibration rotation 
transactions when it is excited. The main mechanism for excitation of the 
vibrational levels of NO'*' is the radiation from the earth at night and the 
earth plus the sun in the daytime. The second mechanism is collisional excita- 
tion. Since the vibrational energy of NO"*" is within 18 cm“^ of the energy of 
N 2 , the cross section for interchange of vibrational energy with N 2 is large. 
The molecule N 2 is infrared inactive and its vibrational temperature should 
exceed the gas kinetic temperature. The reaction rate for collisional excita- 
tion by N 2 was estimated from other resonant and near resonant reactions with 
N 2 . The resulting vibrational temperatures for NO"^ vary between 220° K at 
night at 400 km altitude where collisions are negligible, to 290° K during the 
day at 400 km. Collisional excitations bring the excitation temperature to 
over 360° K at 200 km when T = 2000° K. The average vibrational temperature 
is around 300° K. NO"*” vibration rotation intensities are given in table 14 
for T = 940° K. The column densities needed for the calculations in table 14 
were taken from table 12. However, it should be noted here that the densities 
of N0“^ found by Goldberg and Blumle (ref. 32) for T ~ 1100° K are less than 
the densities of NO"^ given by Hoffman et at. (ref. 30) at T = 940° K for 
altitudes less than 300 km. For lower gas temperatures, the measured density 
of NO'*’ IS much less. 
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TABLE 14.- VIBRATION-ROTATION FUNDAMENTAL OF NO+ 
AT A = 4.27y; 7 = 940° K 


Altitude, 

km 

Intensity, ergs 

-2 -1 -1 
cm ^ sec sr 

Night 

Day 

400 

1.3x10“® 

5.0x10“® 

350 

2.9x10“® 

1.0x10“^ 

300 

1.5x10“® 

3. Ixl0“7 

250 

2.6x10“^ 

1.5x10“® 

200 

4.0x10“® 

8.2x10“® 


The intensities of NO"*” may be higher than the values in table 14 by as 
much as a factor of 200. When a molecule is formed, most of the energy of 
the reaction goes into vibrational excitation of the molecule. The rate of 
formation of NO"*" is about 10^ larger than the rate of excitation of the first 
vibrational level from collisions and radiative excitations. Consequently, 
the rate of spontaneous emission will be larger than that given in table 14 
depending on the fraction of NO"*" molecules formed in an excited state. (The 
rate of formation of NO is not fast enough for this excitation mechanism to 
be important.) 

The intensity of NO”^ will be much stronger for higher temperatures. It 
is thought that the most important reaction for forming NO"*" is 
O"*" + N 2 NO"*" + N, with N 2 '*' + 0 -> NO'*’ + N almost as important. If the abun- 
dance of NO"*" scales as N 2 , then the intensity of the 4.27y vibration-rotation 
fundamental at 350 km should be 50 percent larger at 7 = 1000° K and larger 
by a factor of 50 for 7 = 2000° K. Actually, the density of NO"*" increases 
faster with temperature than the density of N 2 . The reason is that the cross 
section for the reaction O’*" + N 2 NO + N increases with increasing vibra- 
tional temperature of N 2 . Thus at higher gas temperatures, more NO'*' is formed 
per N 2 molecule. 

The rotational line strengths of NO'*' were calculated using the gas tem- 
perature because the rotational energies are coupled much more strongly to 
the translational energies. At 7 = 1000° K, the strongest line is at 87y and 
has an intensity of 1.3x1 O'”® ergs cm“^ sec“^ sr~^ at 350 km. However, at 
7 = 2000° K, the strongest line would be at 60y and would have an intensity 
of about 1.3x10“® ergs cm“^ sec“^ sr~^. 


AN EXAMPLE 


It is necessary to know the infrared fluxes from the upper atmosphere 
when planning infrared astronomical telescopes to be carried on satellites 
such as the space shuttle. It is currently planned that the space shuttle 
infrared telescope will be 1 m in diameter and a typical detector will have a 
1-arc-min field of view. The most interesting quantity is the flux of 
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energy E entering the detector. For a perfectly efficient system, this is 
given by 

E = lAQ. 


where A is the area of the aperture of the telescope and is the solid angle 
of the beam, A 1-m telescope has A = 0.785x10^^ cm^. A beam of 1 arc min^ 
has = 0.846x10“^ sr. Thus a 1-m telescope with a beam of 1 arc min^ will 
measure J x 0.6646x10“^® W if J is in ergs cm“^ sec“^ sr“^. 

The noise due to the photon fluctuations of the radiation can be expressed 
as NEP = X Ep where Ep is the energy per photon = hv . For the above 1-m 

telescope, NEP = /(I . 320x10“^^ J/Xy) W/Hz^/^ where I is in 
ergs cm~^ sec~^ ster”^ and ly is the wavelength in microns. 

A summary of these calculations for a 1-m telescope with a beam of 
one arc minute at 350 km when T = 1000° K on the day side of the earth is 
given in table 15. The NO and NO"*" rotation lines are summed into lOp bands. 

TABLE 15.- RADIATION FROM THE ATMOSPHERE ENTERING A 1-m TELESCOPE WITH A 


1-arc-min FIELD OF VIEW 


A, M 

Watts into the detector 

NEP, W/Hz-1/2 

[01] 63.1 

8.0x10"^5 

5.0x10-10 

[01] 147 

2.0x10"^6 

5.2xlO~10 

[N II] 121.7 

6.0x10-20 

9.9x10-21 

[OIII] 51.71 

6.3x10-20 

1.6x10-20 

O 2 7''. 8 

5.3x10-22 

1.2x10-21 

CO 4.7 

l.OxlO-1^ 

6.5x10-10 

NO 5.3 

2.3xl0“10 

2.9x10-10 

N0+ 4.3 

1.0x10"^^ 

7.0x10-10 

CO + NO + N0+ 60-70 

8.0x10-10 

1.6x10-10 

CO + NO + N0+ 70-80 

1.0x10-1^ 

1.6x10-10 

CO + NO + N0+ 80-90 

1.0x10-11' 

1.5x10-10 


The line strengths will Increase by factors of 2 to 5 for every 50 km that 
one goes down in altitude, and will Increase by factors of 5 to 10 as the 
temperature of the gas is doubled. One must remember that T = 1000° K is 
only an average temperature, and that 1 percent of the time the temperature 
is over 1900° K. In addition, the path length and the column density increase 
as one looks away from the zenith. Except for very large zenith angles 0, 
the column density and the path length increase as sec 0; this is a factor 
of 2 for 0 = 60°. 


SUMMARY OF RESULTS 


Infrared radiation from atoms and molecules in the earth's upper atmo- 
sphere has been calculated for zenith lines-of-sight at altitudes from 200 km to 
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400 km. The infrared spectra at each of five altitudes are plotted in fig- 
ures l(a)-(e). Because of the great uncertainty in extrapolating the abun- 
dances, only spectra for the average condition T = 1000° K have been plotted. 
However, rough estimates can be made of the intensities for T = 2000° K. The 
intensity of the NO"*" 4.3-y band at 300 km is increased by a factor of 25 on 
the sunlit side and by a factor of 200 on the dark side (i.e., collisional 
excitation now dominates over radiative excitation), the NO band at 5.3y is 
increased by a factor of 25, and the CO band at 4.7y is increased by a factor 
of 70 on the sunlit side and by a factor of 150 on the dark side. The energy 
in the NO + NO"^ + CO pure rotation lines is increased by a factor of 50 and 
the maximum emission is in the 50y-60y band Instead of 80y-90y band. 

Because accurate abundances and excitation cross sections were not 
available for many of the cases considered, the predictions of infrared 
intensity are only estimates. For neutral atomic oxygen, the intensities 
for all temperatures are probably good to 50 percent, but for the ions and 
molecules the abundances and column densities are much less certain. Compu- 
tations of the NO and CO abundances have only been done for average tempera- 
tures. At higher temperatures the Increased density will make reactions of 
formation and destruction more likely, and thus make the abundance of NO and 
CO scale differently from that given by pure diffusion (i.e., like N 2 ) . 
Similarly, higher temperatures are caused by increased solar UV flux, which 
in turn causes increased ionization, but higher densities bring about increased 
recombination. For the vibration-rotation bands of NO"*" there is additional 
uncertainty because the band strengths have not been measured in the labora- 
tory. Moreover, the collisional excitation rates for all the molecules are 
quite uncertain. The final uncertainties (95 percent confidence) for 
T = 1000° K might be as high as one order of magnitude and minus one-half 
order of magnitude for the NO"^ vibration- rotation band, plus and minus one- 
half order of magnitude for the NO and CO vibration-rotation bands, and plus 
one-half order of magnitude and minus one order of magnitude for the pure 
rotation lines. The uncertainties at 7 = 2000° K are one-half an order of 
magnitude larger in the molecular lines. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Oct. 23, 1975 
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